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Measurement of prompt J/i/a and beauty hadron production cross sections 
at mid-rapidity in pp collisions at i/s = 7 TeV 



The ALICE experiment at the LHC has studied J/y/ production at mid-rapidity in pp collisions at 
y/s = 7 TeV through its electron pair decay on a data sample corresponding to an integrated lu- 
minosity L; nt = 5.6nb _1 . The fraction of J/y from the decay of long-lived beauty hadrons was 
determined for J/y candidates with transverse momentum p t > 1.3 GeV/c and rapidity \y\ < 0.9. 
The cross section for prompt J/y/ mesons, i.e. directly produced J/y and prompt decays of heavier 
charmonium states such as the y(2S) and Xc resonances, is CJpmmpts/yiPt > 1-3 GeV/c, \y\ < 0.9) 
= 8.3 ±0.8 (stat.) ±1.1 (syst.) 4 (syst.pol.) /ib. The cross section for the production of b-hadrons 
decaying to J/i/a with p t > 1.3 GeV/c and \y\ < 0.9 is Oj/ ¥< _ bB {Pt > 1-3 GeV/c, \y\ < 0.9) = 1.46 
± 0.38 (stat.) ^032 (syst.) /lb. The results are compared to QCD model predictions. The shape of 
the pt and y distributions of b-quarks predicted by perturbative QCD model calculations are used to 
extrapolate the measured cross section to derive the bb pair total cross section and dc/dy at mid- 
rapidity. 
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1 Introduction 

The production of both charmonium mesons and beauty-flavoured hadrons, referred to as b-hadrons or 
he in this paper, in hadronic interactions represents a challenging testing ground for models based on 
Quantum ChromoDynamics (QCD). 

The mechanisms of J/y production operate at the boundary of the perturbative and non-perturbative 
regimes of QCD. At hadron colliders, J/y production was extensively studied at the Tevatron 03-@l and 
RHIC Q. Measurements in the new energy domain of the Large Hadron Collider (LHC) can contribute 
to a deeper understanding of the physics of the hadroproduction processes. The first LHC experimen- 
tal results on the J/i//" transverse momentum (p t ) differential cross sections rtfrj-TTOl are well described 
by various theoretical approaches ITIT1 - PT41 . Among those results, the ALICE Collaboration reported the 
measurement of the rapidity (y) and transverse momentum dependence of inclusive J/y production in 
proton-proton (pp) collisions at y/s = 7 TeV (9J. The inclusive J/y yield is composed of three contribu- 
tions: prompt J/y produced directly in the proton-proton collision, prompt J/y produced indirectly (via 
the decay of heavier charmonium states such as % c and yr(2S)), and non-prompt J/y from the decay of 
b-hadrons. Other LHC experiments have separated the prompt and non-prompt J/y component |[6 K8lfT0l . 
However, at mid-rapidity, only the high-/? t part of the differential dO]/y/dp t distribution was measured 
(p t > 6.5 GeV/c), i.e. a small fraction (few percent) of the /^-integrated cross section. 

The measurement of the production of b-hadrons in pp collisions at the LHC provides a way to test, in a 
new energy domain, calculations of QCD processes based on the factorization approach. In this scheme, 
the cross sections are computed as a convolution of the partem distribution functions of the incoming 
protons, the partonic hard scattering cross sections, and the fragmentation functions. Measurements of 
cross sections for beauty quark production in high-energy hadronic interactions have been done in the 
past at pp colliders at center-of-mass energies from 630 GeV lfT51,[T6l to 1.96 TeV |[2l [T7HT9l and in p- 
nucleus collisions with beam energies from 800 to 920 GeV ||20"1 . The LHC experiments have reported 
measurements of b-hadron production in pp collisions at y/s = 7 TeV by studying either exclusive decays 
of B mesons II2D4231 or semi-inclusive decays of b-hadrons ll6]48l [T0ll24ll25l . At mid-rapidity, the mea- 
surements are available only for p t of the b-hadrons larger than 5 GeV/c, whereas the low p t region 
of the differential b-hadron cross sections, where the bulk of the b-hadrons is produced, has not been 
studied. 

In this paper, the measurement of the fraction of J/y from the decay of b-hadrons in pp collisions at 
■yfs = 7 TeV for J/y in the ranges 1.3 < p t < 10 GeV/c and \y\ < 0.9 is determined. This information 
is combined with the previous inclusive J/y cross section measurement reported by ALICE ||9). Prompt 
J/y and b-hadron cross sections are thus determined at mid-rapidity down to the lowest p t reach at the 
LHC energy. 

2 Experiment and data analysis 

The ALICE experiment |[26l consists of a central barrel, covering the pseudorapidity region | rj | < 0.9, 
and a muon spectrometer with — 4 < r\ < — 2.5 coverage. The results presented in this paper were 
obtained with the central barrel tracking detectors, in particular the Inner Tracking System (ITS) 1126 11271 
and the Time Projection Chamber (TPC) ESI . The ITS, which consists of two innermost Silicon Pixel 
Detector (SPD), two Silicon Drift Detector (SDD), and two outer Silicon Strip Detector (SSD) layers, 
provides up to six space points (hits) for each track. The TPC is a large cylindrical drift detector with 
an active volume that extends over the ranges 85 < r < 247 cm and —250 < z < 250 cm in the radial 
and longitudinal (beam) directions, respectively. The TPC provides up to 159 space points per track and 
charged particle identification via specific energy loss (dE / dx) measurement. 

The event sample, corresponding to 3.5 x 10 8 minimum bias events and an integrated luminosity Lj nt = 
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5.6nb _1 , event selection and track quality cuts used for the measurement of the inclusive J/y production 
at mid-rapidity [9 ] were also adopted in this analysis. In particular, an event with a reconstructed vertex 
position z v was accepted if |z v | < 10 cm. The tracks were required to have a minimum p t of 1 GeV/c, a 
minimum number of 70 TPC space points, a % 2 per space point of the momentum fit lower than 4, and 
to point back to the interaction vertex within 1 cm in the transverse plane. At least one hit in either of the 
two layers of the SPD was required. For tracks passing this selection, the average number of hits in the 
six ITS layers was 4.5-4.7, depending on the data taking period. The electron identification was based 
on the specific energy loss in the TPC: a ±3a inclusion cut around the Bethe-Bloch fit for electrons 
and ±3. 5a (±3a) exclusion cut for pions (protons) were employed ||9]. Finally, electron or positron 
candidates compatible, together with an opposite charge candidate, with being products of y conversions 
(the invariant mass of the pair being smaller than 100 MeV/c 2 ) were removed, in order to reduce the 
combinatorial background. It was verified, using a Monte Carlo simulation, that this procedure does not 
affect the J/y signal. In this analysis, opposite-sign (OS) electron pairs were divided in three "types": 
type "first-first" (FF) corresponds to the case when both the electron and the positron have hits in the 
first pixel layer, type "first-second" (FS) are those pairs where one of them has a hit in the first layer and 
the other does not, while for the type "second-second" (SS) neither of them has a hit in the first layer. 
The candidates of type SS, which correspond to about 10% of the total, were discarded due to the worse 
spatial resolution of the associated decay vertex. 



A detailed description of the track and vertex reconstruction procedures can be found in [29|. The pri- 
mary vertex was determined via an analytic x 2 minimization method in which tracks are approximated as 
straight lines after propagation to their common point of closest approach. The vertex fit was constrained 
in the transverse plane using the information on the position and spread of the luminous region. The 
latter was determined from the distribution of primary vertices reconstructed over the run. Typically, the 
transverse position of the vertex has a resolution that ranges from 40 jltm in low-multiplicity events with 
less than 10 charged particles per unit of rapidity to about 10 /im in events with a multiplicity of about 
40. For each J/y candidate a specific primaiy vertex was also calculated by excluding the J/y decay 
tracks, in order to estimate a systematic uncertainty related to the evaluation of the primary vertex in the 
case of events with non-prompt J/y, as discussed in section[3] The decay vertex of the J/y candidate was 
computed with the same analytic % 2 minimization as for the primary vertex, using the two decay tracks 
only and without the constraint of the luminous region. 

The measurement of the fraction of the J/y yield coming from b-hadron decays, f&, relies on the dis- 
crimination of J/y mesons produced at a distance from the pp collision vertex. The signed projection 
of the J/y flight distance onto its transverse momentum vector, p t , was constructed according to the 
formula 

^ = 1-^/^, (1) 

where L is the vector from the primary vertex to the J/y decay vertex. The variable x, referred to 
as "pseudoproper decay length" in the following, was introduced to separate prompt J/y from those 
produced by the decay of b-hadronfl 



j/y 
Pt 



(2) 



where m } / ¥ is the (world average) J/y mass 113031- 



For events with very low J/y p t , the non-negligible amount of J/y with large opening angle between its 
flight direction and that of the b-hadron impairs the separation ability. Monte Carlo simulation shows 



1 The variable x, which was introduced in [1], mimics a similar variable used for b-hadron lifetime measurements where 
b-hadrons are reconstructed exclusively and therefore the mass and p t of the b-hadron can be used in place of those of the J/y/; 
to get c% = jj = V„L ■ 
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that the detector resolution allows the determination of the fraction of J/y from the decay of b-hadrons 
for events with J/y p t greater than 1.3 GeV/c. 

An unbinned 2-dimensional likelihood fit was used to determine the ratio of the non-prompt to inclu- 
sive J/y production and the ratio of J/y signal candidates (the sum of both prompt and non-prompt 
components) to the total number of candidates, /sig> by maximizing the quantity 

N 

\nL = Y l foF(x,m ft + e -), (3) 

where m e + e - is the invariant mass of the electron pair and N is the total number of candidates in the range 
2.4 < n? e + e - < 4.0 GeV/c 2 . The expression for F(jc,m e + e -) is 

F(x,m e+e -) = / sig -F S i g (x) ■ Afsig(m e + e -) + (1 - fsig) ■ F Bks (x) ■ M Bkg (m s+e -), (4) 

where FsigC*) and F Bkg (x) are Probability Density Functions (PDFs) describing the pseudoproper decay 
length distribution for signal and background candidates, respectively. Msi g (m e + e - ) and M B ^ g (m e + e - ) are 
the PDFs describing the dielectron invariant mass distributions for the signal and background, respec- 
tively. A Crystal Ball function OTTl is used for the former and an exponential function for the latter. The 
signal PDF is given by 

^ig (x) = f B ■ F B (x) + ( 1 - f B ) • Fprompt (x) , (5) 

where F pmmpt (x) and F B (x) are the PDFs for prompt and non-prompt J/y, respectively, and f B is the 
fraction of reconstructed non-prompt J/y, 

jB = 77 -T7 . (O) 

■ ,v J/V/^h B "•""promptJ/v/ 

which can differ (see below) from f B due to different acceptance and reconstruction efficiency of prompt 
and non-prompt J/y. The distribution of non-prompt J/y is the convolution of the x distribution of J/y 
from b-hadron events, % B (x), and the experimental resolution on x, R type (x), which depends on the type 
of candidate (FF or FS), 

F B (x) = XB (X 1 ) ® Rtype ~ x) . (7) 

Promptly produced J/y mesons decay at the primary vertex, and their pseudoproper decay length distri- 
bution is thus simply described by R type (x): 

Fprompt (x) = S(x')<g> R type (x 1 ~x)= R type (x) . (8) 

The resolution function is described by the sum of two Gaussians and a power law function reflected 
about x = and was determined, as a function of the p t of the J/y, with a Monte Carlo simulation 
study. In this simulation, which utilizes GEANT3 02l and incorporates a detailed description of the 
detector material, geometry, and response, prompt J/y were generated with a p t distribution extrapo- 
lated from CDF measurements [T| and a y distribution parameterization taken from Color Evaporation 
Model (CEM) calculations f33l . These J/y were individually injected into proton-proton collisions sim- 
ulated using the PYTHIA 6.4.21 event generator 0411351 . and reconstructed as for J/y candidates in 
data. A data-driven method (discussed in section [3} was also developed and used to estimate the sys- 
tematic uncertainty related to this procedure. The Monte Carlo x distribution of J/y from the decay 
of b-hadrons produced in proton-proton collisions simulated using the PYTHIA 6.4.21 event genera- 
tor 1134 11351 with Perugia-0 tuning Il36l was taken as the template for the x distribution of b-hadron events 
in data, % B {x). A second template, used to estimate the systematic uncertainty, was obtained by de- 
caying the simulated b-hadrons using the EvtGen package 11371 . and describing the final state radiation 
("internal" bremsstrahlung) using PHOTOS 
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Fig. 1: Invariant mass (left panel) and pseudoproper decay length (right panel) distributions of opposite sign 
electron pairs for \yj/y\ < 0.9 and > 1-3 GeV/c w hh superimposed projections of the maximum likelihood fit. 
The latter distribution is limited to the i/y/ candidates under the mass peak, i.e. for 2.92 < m e + e - < 3.16 GeV/c 2 , 
for display purposes only. The % 2 values of these projections are reported for both distributions. 

For the background x distribution, F^ g {x), the functional form employed by CDF O] was used, 



^BkgW =(1 -/- 
+ 



f— fsymjRtype {%) 



sym 



MA, 



21 



sym 



(9) 



where 6(x) is the step function, /+, /_ and / sym are the fractions of three components with positive, neg- 
ative and symmetric decay length exponential distributions, respectively. The effective parameters A + , 
A_ and A sym , and optionally also the corresponding fractions, were determined, prior to the likelihood fit 
maximization, with a fit to the x distribution in the sidebands of the dielectron invariant mass distribu- 
tion, defined as the regions 1.8-2.6 and 3.2-5.0 GeV/c 2 . The introduction of these components is needed 
because the background consists also of random combinations of electrons from semi-leptonic decays 
of charm and beauty hadrons, which tend to produce positive x values, as well as of other secondary or 
mis-reconstructed tracks which contribute both to positive and negative x values. The first term in eq. [9] 
proportional to R type {x), describes the residual combinatorics of primary particles. 

In figure [Q the distributions of the invariant mass and the pseudoproper decay length, the latter restricted 
to candidates with 2.92 < m e + e - < 3.16 GeV/c 2 , for opposite-sign electron pairs with p t > 1.3 GeV/c 
are shown with superimposed projections of the maximum likelihood fit result. 

The value of the fit parameter /g provides the fraction of non-prompt J/y which were reconstructed. 
In principle prompt and non-prompt J/\j/ can have different acceptance times efficiency (A x e) values. 
This can happen because of two effects: (i) the Axe depends on the p t of the J/\j/ and prompt and 
non-prompt J/i/a have different p t distributions within the considered p t range; (ii) at a given p t , prompt 
and non-prompt J/\j/ can have different polarization and, therefore, a different acceptance. The fraction 
of non-prompt J/y/\ corrected for these effects, was obtained as 



,/b 



i | l-f B (Axe) B V 1 

/j3 X *V prompt/ 



(10) 



where (A x e) B and (A x £)p rom pt are the average acceptance times efficiency values, in the considered 
p t range and for the assumed polarization state, of non-prompt and prompt J/\j/, respectively. The accep- 
tance times efficiency (A x e) varies very smoothly with p t and, for unpolarized J/\j/ in the p t range from 
1.3 to 10 GeV/c, has a minimum of 8% at 2 GeV/c and a broad maximum of 12% at 7 GeV/c |9|. As a 
consequence, the (A x e) values of prompt and non-prompt J/y differ by about 3% only in this integrated 
Pt range. 
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The central values of the resulting cross sections are quoted assuming both prompt and non-prompt J/y 
to be unpolarized and the variations due to different assumptions are estimated as a separate systematic 
uncertainty. The polarization of J/y from b-hadron decays is expected to be much smaller than for 
prompt J/y due to the averaging effect caused by the admixture of various exclusive B — > J/y+X decay 
channels. In fact, the sizeable polarization, which is observed when the polarization axis refers to the 
B-meson direction 1401 . is strongly smeared when calculated with respect to the direction of the daughter 
J/y EL as indeed observed by CDF 0. Therefore, these variations will be calculated in the two cases 
of prompt J/y with fully transverse (A = 1) or longitudinal (A = — 1) polarization, in the Collins-Soper 
(CS) and helicity (HE) reference frames^, the non-prompt component being left unpolarized. 

Despite the small J/y candidate yield, amounting to about 400 counts, the data sample could be divided 
into four p t bins (1.3-3, 3-5, 5-7 and 7-10 GeV/c), and the fraction /b was evaluated in each of them 
with the same technique. At low p t the statistics is higher, but the resolution is worse and the signal over 
background, S/B, is smaller (i.e. /sig is smaller). At high p t the statistics is smaller, but the resolution 
improves and the background becomes negligible. In figure |2]the distributions of the invariant mass and 
of the pseudoproper decay length are shown in different p t bins with superimposed results of the fits. 

3 Systematic uncertainties 

The different contributions to the systematic uncertainties affecting the measurement of the fraction of 
J/y from the decay of b-hadrons are discussed in the following, referring to the integrated p t range, and 
summarized in table Q] 

- Resolution function. The resolution function was determined from a Monte Carlo simulation, as 
discussed above. The fits were repeated by artificially modifying the resolution function, according 
to the formula 



where 8 is a constant representing the desired relative variation of the RMS of the resolution 
function. Studies on track distance of closest approach to the primary interaction vertex in the 
bending plane (do) show that the p t dependence of the do resolution as measured in the data is 
reproduced within about 10% by the Monte Carlo simulation ll29l . but with a systematically worse 
resolution in data. For the x variable a similar direct comparison to data is not straightforward, 
however, the residual discrepancy is not expected to be larger than that observed for do- 

The variations of /b obtained in the likelihood fit results by varying 8 from —5% to +10% are 
+8% and -15%, respectively, and they were assumed as the systematic uncertainty due to this 
contribution. 

An alternative, data-driven, approach was also considered. The x distribution of the signal, com- 
posed of prompt and non-prompt J/y, was obtained by subtracting the x distribution of the back- 
ground, measured in the sidebands of the invariant mass distribution. This distribution is then fitted 
by fixing the ratio of prompt to non-prompt J/y to that obtained from the likelihood fit and leaving 
free the parameters of the resolution function. The RMS of the fitted resolution function is found 
to be 8% larger than the one determined using the Monte Carlo simulation, hence within the range 
of variation assumed for 8. 

- Pseudoproper decay length distribution of background. The shape of the combinatorial back- 
ground was determined from a fit to the x distribution of candidates in the sidebands of the invariant 
mass distribution. By varying the fit parameters within their errors an envelope of distributions was 

2 The polar angle distribution of the J/yr decay leptons is given by cW/dcos 6 = 1 + A cos 2 6. 
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Fig. 2: Invariant mass (left panels) and pseudoproper decay length (right panels) distributions in different p t bins 
with superimposed projections of the maximum likelihood fit. The % 2 values of these projections are also reported 
for all distributions. 

obtained, whose extremes were used in the likelihood fit in place of the most probable distribution. 
The variations in the result of the fit were determined and adopted as systematic uncertainties. 
Also, it was verified that the x distribution obtained for like-sign (LS) candidates, with invariant 
mass in the range from 2.92 to 3.16 GeV/c 2 complementary to the sidebands, is best fitted by a 
distribution which falls within the envelope of the OS distributions. Finally, the likelihood fit was 
repeated by relaxing, one at a time, the parameters of the functional form (eq. [9} and it was found 
that the values of /b were within the estimated uncertainties. The estimated systematic uncertainty 
is 6%. 

- Pseudoproper decay length distribution of b-hadrons. The fits were also done using as template 
for the x distribution of b-hadrons, Xb (x) , that obtained by the EvtGen package 11371 , and describing 
the final state radiation using PHOTOS [38,39]. The central values of the fits differ by a few percent 
at most and the resulting systematic uncertainty is 3%. 

- Invariant mass distributions. The likelihood method was used in this analysis to fit simulta- 
neously the invariant mass distribution, which is sensitive to the ratio of signal to all candidates 
(/sig)> an d the x distribution, which determines the ratio of non-prompt to signal candidates (/b)- 
The statistical uncertainties on these quantities were therefore evaluated together, including the 
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effects of correlations. However, the choice of the function describing the invariant mass distri- 
bution, as well as the procedure, can introduce systematic uncertainties in the evaluation of /b- 
Different approaches were therefore considered: (i) the functional form describing the background 
was changed into an exponential plus a constant and the fit repeated; (ii) the background was de- 
scribed using the LS distribution and the signal was obtained by subtracting the LS from the OS 
distributions. The signal and the background shapes were determined with % 2 minimizations. Both 
functional forms, exponential and exponential plus a constant, were considered for the background. 
The likelihood fit was then performed again to determine fg (and /si g ); (Hi) the same procedure 
as in ( ii) was used, but additionally /sig was estimated a priori using a bin counting method @ 
instead of the integrals of the best fit functions. The maximum likelihood fit was performed with 
/sig fixed to this new value; (iv) and (v) the same procedures as in (ii) and (Hi) were used but with 
the background described by a track rotation (TR) method ||9] . 

Half of the difference between the maximum and minimum /b values obtained with the different 
methods was assumed as systematic uncertainty. It amounts to about 6%. 

- Primary vertex. The effect of excluding the decay tracks of the J/\j/ candidate in the computation 
of the primary vertex was studied with the Monte Carlo simulation: on the one hand, for the prompt 
J/y, the x resolution function is degraded, due to the fact that two prompt tracks are not used in the 
computation of the vertex, which is thus determined with less accuracy. The effect on the resolution 
is p t dependent, with the RMS of the x distribution of prompt J/i//" increasing by 15% at low p t and 
by 7% at high p t . On the other hand, for non-prompt J/y a bias on the x determination should be 
reduced. The bias consists in an average shift of the primary vertex towards the secondary decay 
vertex of the b-hadrons, which is reflected in a shift of the mean of the x distribution by about 4 jUm 
for the /^-integrated distribution. However, the shift is p t and "type" dependent. In some cases 
the bias is observed in the opposite direction and is enhanced by removing the decay tracks of the 
candidate. This can happen since b-quarks are always produced in pairs. If a charged track from 
the fragmentation of the second b-quark also enters the acceptance, it can pull the primary vertex 
position towards the opposite direction. In the end, therefore, the primaiy vertex was computed 
without removing the decay tracks of the candidates. To estimate the systematic uncertainty, the 
analysis was repeated by either ( i) removing the decay tracks in the computation of the primary 
vertex and using the corresponding worse resolution function in the fit or (ii) keeping those tracks 
and introducing an ad hoc shift in the distribution of the Xb (*)> equal to that observed in the Monte 
Carlo simulation for non-prompt J/\j/. The contribution to the systematic uncertainty is about 5%. 

- MC p t spectrum. The ratio ,}** s e)b in eq.Q! was computed using MC simulations: prompt J/w 

[AXL) p ro mpt 

were generated with the p t distribution extrapolated from CDF measurements [Q] and the y distri- 
bution parameterized from CEM D~3~l : b-hadrons were generated using the PYTHIA 6.4.21 041351 
event generator with Perugia-0 tuning [36 ]. By varying the average p t of the J/\j/ distributions 
within a factor 2, a 1.5% variation in the acceptance was obtained both for prompt and non-prompt 
J/\j/. Such a small value is a consequence of the weak p t dependence of the acceptance. For the 
measurement integrated over p t (p t > 1.3 GeV/c), the Axe values of prompt and non-prompt J/\j/ 
differ by about 3% only. The uncertainty due to Monte Carlo p t distributions is thus estimated to 
be 1%. When estimating /b in Pt bins, this uncertainty is negligible. 

- Polarization. The variations of /b obtained assuming different polarization scenarios for the 
prompt component only were evaluated, as discussed in section |2l and are reported in table Q] 
The maximum variations are quoted as separate errors. 

The study of systematic uncertainties was repeated as a function of p t . In table Q] the results are sum- 
marized for the integrated p t range (p t > 1.3 GeV/c) and for the lowest (1.3-3 GeV/c) and highest 
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Table 1: Systematic uncertainties (in percent) on the measurement of the fraction of J/y/ from the decay of b- 
hadrons, /g. The variations of /g are a ls° reported, with respect to the case of both prompt and non-prompt J/y/ 
unpolarized, when assuming the prompt component with given polarization. 



Source 


Systematic uncertainty (%) 




Pt integrated 


lowest p t bin 


highest p t bin 


Resolution function 


+8, -15 


+15, -25 


+2, -3 


x distribution of background 


±6 


±13 


±1 


x distribution of b-hadrons 


±3 


±3 


±2 


m e + e - distributions 


±6 


±11 


±4 


Primary vertex 


+4, -5 


±4 


+4, -8 


MC p t spectrum 


±1 








Total 


+12,-18 


+23, -30 


+6,-9 


Polarization (prompt J/y/") 








CS (A = -l) 


+13 


+22 


+5 


CS (A = +1) 


-10 


-19 


-3 


HE (A = -1) 


+17 


+19 


+11 


HE (A = +1) 


-14 


-16 


-8 



(7-10 GeV/c) pt bins. All systematic uncertainties increase with decreasing p t , except the one related to 
the primary vertex measurement. 

4 Results 

4.1 Fraction of J/y/ from the decay of b-hadrons 

The fraction of J/y/ from the decay of b-hadrons in the experimentally accessible kinematic range, p t > 
1.3 GeV/c and \y\ < 0.9, which is referred to as "measured region" in the following, is 

h = 0.149±0.037(stat.)+°;^ ^0^21 (t^-i) (syst.pol.). 

The fractions measured in the p t bins are reported in tableland shown in figure [3] In the figure, the data 
symbols are placed at the average value of the p t distribution of each bin. The average was computed 
using the above mentioned Monte Carlo distributions: the one based on the CDF extrapolation ll33l and 
that using PYTHIA (34l[35l with Perugia-0 tuning 11361 for prompt and non-prompt J/y/, respectively, 
weighted by the measured /b- In figure [3] the results of the ATLAS [8] and CMS ifTOl experiments 
measured at mid-rapidity for the same colliding system are also shown. The ALICE results extend the 
mid-rapidity measurements down to low p t . 

4.2 Prompt J/\j/ production 

By combining the measurement of the inclusive J/y/ cross section, which was determined as described 
in (9], and the /b value, the prompt J/y/ cross section was obtained: 

°promptJ/v/ = (1 — /b) • Oj/y (H) 

The numerical values of the inclusive J/y/ cross section in the p t ranges used for this analysis are sum- 
marized in table [2] In the measured region the integrated cross section is CJ prompt j/y(|;y| < 0.9, p t > 

1.3GeV/c) = 8.3 ±0.8(stat.) ± l.l(syst.)+j ^j^ 1 ^ /xb. The systematic uncertainties related to the un- 
known polarization are quoted for the reference frame where they are the largest. 

The differential distribution is shown as a function of p t in figure g] and d< W J /r is p i ott ed in 

figure [5] The numerical values are summarized in table |2] In figure [4] the statistical and all systematic 
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Fig. 3: The fraction of J/y/ from the decay of b-hadrons as a function of p t of J/y/ compared with results from 
ATLAS HD and CMS QUI in pp collisions at y/s =7 TeV. 

errors are added in quadrature for better visibility, while in figure [5] the error bar shows the quadratic 
sum of statistical and systematic errors, except for the 3.5% systematic uncertainty on luminosity and 
the 1% on the branching ratio (BR), which are added in quadrature and shown as box. The results shown 
in figures |4] and [5] assume unpolarized J/y production. Systematic uncertainties due to the unknown 
J/y polarization are not shown. Results by the CMS EES, LHCb Q and ATLAS HI Collaborations 
are shown for comparison. Also for these data the uncertainties due to luminosity and to the BR are 
shown separately (boxes) in figure [5] while the error bars represent the statistical and the other sources 
of systematic uncertainties added in quadrature. 



The ALICE 



prompt S/y 

dydpt 



measurement at mid-rapidity (left panel of figure 5J is complementary to the data 
of CMS, available for \y\ < 0.9 and p t > 8 GeV/c, and ATLAS, which covers the region |y| < 0.75 and 
Pt> 1 GeV/c. In the right panel of figure 0] the ALICE results are compared to next-to-leading order 
(NLO) non-relativistic QCD (NRQCD) theoretical calculations by M. Butenschon and B.A. Kniehl fT2ll 
and Y.-Q. Ma et al. ifTBl . Both calculations include color-singlet (CS), color-octet (CO), and heavier 
charmonium feed-down contributions. For one of the two models (M. Butenschon and B.A. Kniehl) 
the partial results with only the CS contribution are also shown. The comparison suggests that the CO 
processes are indispensable to describe the data also at low p t . The results are also compared to the 
model of V.A. Saleev et al. Ifl4ll . which includes the contribution of partonic sub-processes involving 
t-channel parton exchanges and provides a prediction down to p t = 0. 

da™ 



The ALICE result for 



J prompt 



dv 



(figure |5]), which equals 



da. 



prompt J/ v 



= 5.89±0.60(stat.)+^(syst.)+«;[»(extr.)!; : ^^:i\ )M b, 



was obtained by subtracting from the inclusive J/i/a cross section measured for p t > that of J/\j/ coming 
from b-hadron decays. The latter was determined, as discussed in the next section, by extrapolating the 
cross section from the measured region down to p t > using an implementation of pQCD calculations 
at fixed order with next-to leading-log resummation (FONLL) BP . The extrapolation uncertainty is 
negligible with respect to the other systematic uncertainties. In figure [5] the CMS and LHCb results for 
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Table 2: The fraction of J/y/ from the decay of b-hadrons and cross sections. Some of the contributions to the 
systematic uncertainty do not depend on p t , thus affecting only the overall normalization, and they are separately 
quoted (correl.). The contributions which depend on p t , even when they are correlated bin by bin, were included 
among the non-correlated systematic errors. The values of (p t ) were computed using Monte Carlo distributions 
(see text for details). 
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3.07 


CTj/v/^hB (bj/vl < 0-9) (Mb) 
1.46 ±0.38 +0.26,-0.32 











Pt >o 
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1.77 ±0.46 


+0.32, -0.39 
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d<7 b6M'l| v |<0.9 (^ b > 

43±11 +9,-10 +0.6,-1.5 



°bb (^b) 

282 ±74 +58,-68 +8,-7 

the rapidity bins where the p t coverage extends down to zero were selected. For CMS, the value for 
1.6 < |y| < 2.4 was obtained by integrating the published d 2 a prompt j/y/dptdy data The ALICE data 
point at mid-rapidity complements the other LHC measurements of prompt J/i//" production cross section 
as a function of rapidity. It is worth noting that the uncertainties of the data sets of the three experiments 
are uncorrelated, except for that (negligible) of the BR, while within the same experiment most of the 
systematic uncertainties are correlated. The prediction of the model by V.A. Saleev et al. Ifl4l at mid- 
rapidity provides dgp ™ plJ/y = 7.8+4*5 /lb, which, within the large band of theoretical uncertainties, is in 
agreement with our measurement. 

4.3 Beauty hadron production 

The cross section of J/i/ from b-hadrons decay was obtained as Oj/ v< _ hB = /b ■ Oj/y in the measured 
region it is 

Oj/^CPt > 1.3GeV/c, |y| < 0.9) = 1.46±0.38(stat.)^(syst.)/xb. 

This measurement can be compared to theoretical calculations based on the factorization approach. In 
particular, the prediction of the FONLL lETTI . which describes well the beauty production at Tevatron en- 
ergy, provides Ii42l 1 .33^ 48 A 4 b, in good agreement with the measurement. For this calculation CTEQ6.6 
partem distribution functions P3l were used and the theoretical uncertainty was obtained by varying 
the factorization and renormalization scales, jUp and /Ir, independently in the ranges 0.5 < pip/nit < 2, 
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pig _ 4 . dg p^pt Vv as a f unct i on G f pt compared to results from ATLAS fl8] and CMS flU at mid-rapidity (left panel) 
and to theoretical calculations [ 12- 14 1 (right panel). The error bars represent the quadratic sum of the statistical 
and systematic uncertainties. 

0.5 < jUR/m t < 2, with the constraint 0.5 < jUf/Mr < 2, where m t = J p 2 +m^. The beauty quark mass 
was varied within 4.5 < m\, < 5.0 GeV/c 2 . 

The same FONLL calculations were used to extrapolate the cross section of non-prompt J/y down to 
p t equal to zero. The extrapolation factor, which is equal to 1.212^QQ3g, was computed as the ratio of 

the cross section for p J / v > and \yj/w\ < 0.9 to that in the measured region {p\ ,w > 1.3 GeV/c and 
bj/yl < 0.9). Using the PYTHIA event generator with Perugia-0 tuning instead of FONLL provides 
an extrapolation factor of 1.156. The measured cross section corresponds thus to about 80% of the 
/^-integrated cross section at mid-rapidity. Dividing by the rapidity range Ay = 1.8 one obtains 

dCTj /V^h B _ , \+0.18 ^+0.01 



d.v 



0.98 ±0.26 (stat.)™ (syst.)l^ (extr.) jub. 



In figure |6]this measurement is plotted together with the LHCb [7 ] and CMS [Q data at forward rapidity. 
For CMS the values for 1.2 < \y\ < 1.6 and 1.6 < \y\ < 2.4 were obtained by integrating the published 
d 2 Oj/ ¥ ^ hB /dp t dy data J6j; the value for 1.2 < \y\ < 1.6 was also extrapolated from />J™ n = 2.0 GeV/c to 
p t = 0, with the approach based on the FONLL calculations as previously described. The extrapolation 
uncertainties are shown in figure|6]as the slashed areas. The central FONLL prediction and its uncertainty 
band are also shown. A good agreement between data and theory is observed. 



A similar procedure was used to derive the bb quark-pair production cross section 

dcTbb _ d( C 0ry x <rj/ y <-hB(A /y > l-3GeV/c, |y J/y | < 0.9) 
d ^ ^ K <fi£*M N > 1.3GeV/c, \y J/¥ \ < 0.9) 



'J/vr<-h B ' 

where the average branching fraction of inclusive b-hadron decays to J/y measured at LEP I rTTl FT oll . 
BR(hb — > J/y + X) = (1.16 ±0.10)%, was used in the computation of Oj/^ hB - The extrapolation with 
the FONLL calculations provides 

^=43±ll(stat.)t? (syst.)+?j(extr.) fib. 
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Fig. 5: — p'°mp tJ /v as a f unC ];i on G f y The error bars represent the quadratic sum of the statistical and systematic 
errors, while the systematic uncertainties on luminosity and branching ratio are shown as boxes around the data 
points. The symbols are plotted at the center of each bin. The CMS value was obtained by integrating the published 
d 2 <7 p rompt j/vz/dftdy data measured for 1 .6 < \y\ < 2.4 |5J. The results obtained by LHCb f7J and CMS are reflected 
with respect to y — (open symbols). 



Using the PYTHIA event generator with Pemgia-0 tuning (with the EvtGen package to describe the par- 
ticle decays) instead of FONLL results in a central value of 40.4 (40.9) jltb. A compilation of measure- 
ments of da b g/dv at mid-rapidity is plotted in figure UJ as a function of yfs, with superimposed FONLL 
predictions. 

Finally, the total bb cross section was obtained as 

< kk,v^ Oj/^ha^ > 1.3GeV/c, \y iN \ < 0.9) 
a( PP -+ bb +x) = a^-L 2 , g7?(hb ^ J/y + x ; , d3) 

where is the ratio between the yield of J/\j/ mesons (from the decay of b-hadrons) in the full phase 
space and the yield in the measured region \y ]/¥ \ < 0.9 and p\ lw > 1.3 GeV/c. The FONLL calcula- 
tions provide = 4.49+qJq, which produces a(pp ->■ bb+X) = 282±74(stat.)+gg(syst.)+7(extr.) /lb. 
The extrapolation factor was also estimated using PYTHIA with Perugia-0 tuning and found to be 
a 4x TUlA = 4-20. This measurement is in good agreement with those of the LHCb experiment, namely 
288 ±4(stat.) ±48(syst.) ^b and 284 ± 20(stat.) ±49(syst.) /xb, which were based on the measured 
cross sections determined in the forward rapidity range from b-hadron decays into J/t/fX and D°/j.vX, 
respectively 11711241. 



5 Summary 

Results on the production cross section of prompt J/i/a and J/i/a from the decay of b-hadrons at mid- 
rapidity in pp collisions at y/s = 7 TeV have been presented. The measured cross sections have been 
compared to theoretical predictions based on QCD and results from other experiments. Prompt J/\j/ 
production is well described by NLO NRQCD models that include color-octet processes. The cross 
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Fig. 6: — VyftomB as a function of y. The error bars represent the quadratic sum of the statistical and systematic 
errors, while the systematic uncertainties on luminosity and branching ratio are shown as boxes. The systematic 
uncertainties on the extrapolation to p t = are indicated by the slashed areas. The CMS values were obtained 
by integrating the published d 2 Ou V f mmB /dp t dy data measured for 1.2 < \y\ < 1.6 and 1.6 < \y\ < 2.4 J6). The 
results obtained in the forward region by LHCb ||7] are reflected with respect to y = (open symbols). The FONLL 
calculation 14111421 (and its uncertainty) is represented by solid (dashed) lines. 

section of J/\j/ from b-hadron decays is in good agreement with the FONLL prediction, based on per- 
turbative QCD. The ALICE results at mid-rapidity, covering a lower p t region down to p t = 1.3 GeV/c, 
are complementary to those of the ATLAS and CMS experiments, which are available for J/y p t above 
6.5 GeV/c. Using the shape of the p t and y distributions of b-quarks predicted by FONLL calculations, 
the mid-rapidity da /Ay and the total production cross section of bb pairs were determined. 
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